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In this paper we consider dissipative hydrodynamic equations for systems with continuous broken 
symmetries. We first present the case of superfluidity, in which the symmetry U(l) is broken 
and then generalize to the chiral symmetry SU(2)l X SU(2)r. The corresponding new transport 
coefficients are introduced. 



^ ; INTRODUCTION 

o ; 

Relativistic hydrodynamics has often been used as a starting point in the field of heavy ion collisions in order 
to reproduce single-particle spectra or multi-particle correlations (see |jj for example). However, the hot and dense 
hadronic phase created in such experiments is mainly constituted of pions and, until recently [3], no care has been 
taken on the fact that the breaking of the chiral symmetry can affect the theory of hydrodynamics itself 0, 0] . Among 
all the applications of such a modification of the theory, two main directions emerge : first, it is highly wished to have 
a quantitative estimation of this effect when reproducing the experimental spectra j-j and, second, the aim of this 
paper, it is important to complete the whole theory, that is, to include dissipation. 

In hydrodynamic regime, relevant variables are those whose variations in space and time are slow and relaxation 
time becomes infinite in the long wavelength, small frequency limit. Such variables are of two kinds : densities of 
conserved quantities and Goldstone modes associated to continuous broken symmetry [||. The prototype of such a 
£^ , theory which includes both types of hydrodynamic variables is the well-established theory of superfluidity associated 
t^J- ' to the breaking of U(l) and developed many years ago by Landau Q in the non- relativistic case (see also 0, 0| in 
the relativistic domain). The example of the superfluidity will be used in the first part as a guideline in order to 
show explicitly how to include dissipation in hydrodynamic equations with a broken continuous symmetry. Then, in a 
second part, we will generalize to the case of SU(2)l x SU(2)r : starting from the equations of ideal hydrodynamics, 
we will introduce explicitely new transport coefficients associated to the chiral charge densities and Goldstone modes. 



SUPERFLUIDITY 

Let us focus on our first example of hydrodynamics with a continuous broken symmetry : the superfluidity. As 
mentioned in the introduction, the spontaneous breaking of a continuous symmetry implies that the Goldstone modes 
are hydrodynamic variables, i.e. relax very slowly to equilibrium in the long wavelength, small frequency limit. 
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Concerning the superfluid, one therefore has to introduce 0, which is the phase of the condensate which breaks the 
U(l) symmetry associated to the particle number, in our theory. More precisely, the equations for ideal (that is 
without dissipation) relativistic superfluid can be written as |To| : 

d^nou" -V 2 d^4>) = (1) 

g^u = q ( 2 ) 

u^cp = -mo (3) 

where T^ v is the energy-momentum tensor, V 2 is the superfluid density in the non-relativistic limit and /io = 7/1 is 
the chemical potential (7 is the Lorentz factor). We have to notice at this stage that the third equation is equivalent 
to the one contained in the original two-fluid model of Landau |jj when making the correspondence : 

/* l^Landau T" „ ^n-^s 

where v n is the spatial part of the four- velocity = 7(1, v n ) and v s — ^ V(/> is the superfluid velocity. 

We can see that the presence of the new hydrodynamic variable (f> manifests itself not only in the third equation 
but also in the first one (conservation of particle number) and in the second one since T^ v is equal to : 

T» v = (e + p)u>*u" - pg»" + V 2 <j)d v <f> (4) 
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Since we are dealing with an ideal fluid the entropy conservation must be satisfied. Actually it is already contained 
in the second above equation when projected along the direction of u v : u l/ d ll T^ 1 ' = =>■ d^{s a u^) = 0. 

Now, we can go further and add dissipative terms in the hydrodynamic equations. The method, described for 
example in 0, consists in introducing fluxes such as : 

d^nou" - V 2 d"<p + u") = (5) 
((e + P )u^u u - pg» v + V 2 d»4>d v 4> + t"") = (6) 
ufdrf = -Mo - <t>' (?) 

Moreover, our choice for the hydrodynamic velocity u M imposes that u^t^ v = and u^v^ — (the Eckart choice 
leads to other constraints, see jll| for example). 

With these equations, we are now able to derive explicitely the new equation for the entropy. Writing u u d^T^ v = 0, 
we obtain : 

d„(s u» - £V) = -v»dM + ^d,(V 2 d^) + T -^d,u„ (8) 
Jo 1 J o J o 

This equation has the generic form : d^S^ = a where is the total four-flow entropy and a the entropy production. 
a is a bilinear form between fluxes and thermodynamic forces. In the hydrodynamic regime, we are by definition 
near the global equilibrium so that we can express linearly the relation between these fluxes and these forces. The 
coefficients of proportionality are called the transport coefficients. Their physical meaning is thus to characterize the 
magnitude of the response of the system (flows) to a certain disturbance (thermodynamic forces). 

Moreover, because of thermodynamics, a must be positive. This constraint when combined with the Onsager 
principle leads to : 

v» = -uf>u u )d v (£) (9) 

Jo 



dd x (V 2 d x ^ + (C2-lv)dxu x 



+ Tj [{g" x - u"u x )d x u v + {g Xv - u x u v )d x u^)] (10) 
(f>' Q =Cid^ + CMV 2 d^) (11) 



where n is proportional to the thermal conductivity and rj, d, C2, C3 are t ne shear and bulk viscosities (notations are 
those defined in 0). The positivity of a implies that 77, C2, C3 are positive and £ 2 < (2^3 • The sign of £1 has to be 
determined by physical considerations only (here, it can be checked by comparison with [3 that £i is positive as well) . 

CHIRAL DYNAMICS 
Ideal fluid 

In this section, we arc going to recall the main results of hydrodynamics with chiral SU(2) symmetry spontaneously 
broken |2|. The hydrodynamic degrees of freedom in a chiral fluid are the densities of the conserved quantities, namely 
entropy density s, momentum density T , baryonic number density n, left and right-handed charge densities written 
as SU (2)-matrices pl = p l L Ti/2 and p R = p z R Ti/2 and finally the variables associated to the Golstone modes. For chiral 
symmetry, these modes are actually the pions. Then, following 0, we can write the new hydrodynamic variables as 
a 5L/(2)-matrix £ = e lT '' K ^ 7 ' and, by analogy with the superfluid, denote X as "phases". The energy density T 00 is a 
function of all these variables and of the first partial derivatives of £ (assuming £ varies slowly). Following again 
we thus write the energy density T 00 as : 

T 00 = e (s,n,T 0j ) + e 1 

where eo is the "normal fluid" part and ei contains all the non trivial terms, at lowest order, compatible with chiral 
symmetry : 

ex = f -f{5 l0 - 1 - ^ 2 ttV)trasa 3 -St + 2f2n 1 2 2 M p L - Hp R Y?f 
4 l-v 2 v 2 j 2 f 2 (l-v 2 v 2 ) 

1 Mpl + SpRSt) 2 - i v * ? y tr( P£ EG> fc I]t + ^d k Zp R ) 



l 2 fv ' 7 2 (1-«2 V 2) 
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where = f s /ft is the pion velocity and f s ,ft and f v are functions to be determined by thermodynamics of the 
underlying fundamental theory, namely QCD. 

We have fourteen hydrodynamics variables. Therefore we must have fourteen hydrodynamics equations. With the 
hamiltonian density above, it is possible to show 2] that the equations can be written in a covariant way only if we 
make some combinations of the initial variables. Finally the result is : 

d^nou^) = (12) 



dyT^ = 



a,K) + -[Aa] = o 



(13) 
(14) 



^ M ((/ t 2 - f s )u»A + / s 2 £d"Et) + [A, a] = (15) 

with T» v = (e+p)u^u u - pg» v + ^tr^Ed^Et + d y Hd^) and A = u^Ed^E*. 

The nine equations which are specific to chiral dynamics {d^J^ R = and the first order equation for E) are 
actually contained in the first order equation 1)14(1 for a = (pi + Y.pnY^)/^ and in the second order equation 1(15(1 for 
E constructed from the combination (pl — E/9rE^). 

Once again, starting from equation 1(13(1 and using all the other (hydrodynamic and thermodynamic) equations, it 
is possible to deduce the entropy conservation : ^(sou^ 1 ) = as it should be for a perfect fluid. At this stage, we can 
already make some comparisons with the supcrfluid. First of all if E was associated with the breaking of U(l) type 

symmetry, namely E = e 1 ^ , then zEc^E^ would be equal to d M </>. This means that iEVE^ plays the role of superfluid 

f 2 

velocity. Then, by direct comparison between the energy-momentum tensors (for £7(1) and SU(2)), we see that 
can be interpreted as " superfluid density" . 



Dissipation 

We are now in position to treat the dissipation case. As for the superfluid, we add some dissipative flux densities to 
conservation equations : for baryonic number, j£ R = j£ ^j.jTi/2 for chiral charges and t^ v for energy- momentum 
tensor; and we add also a dissipative term for the other non-conserved variables : |£q = |Eq fi/2. 

Then, performing the same combinations as in the previous section, we can finally write after some calculations : 

M (nou" + v") = (16) 
d„ f(e + |))// -pgf*" + ^tr(^E9 l/ E t +8^0^) + t»A = (17) 

d^) + \[A^\a] = -d„? L - E^Et (18) 

id, ((/ t 2 - /?)«M<°> + / 2 E9^Et) + [A<-°\a] = 2d,f L - 2£^j££t (19) 

with A = A^ + |Eg which is the dissipative first order equation for E. A^ represents the non-dissipative part of A. 
With these new equations, it is a simple task to express the equation for the entropy (note that, since the expressions 
for e and p must be given by the same formula as in the non-dissipative case, e and p depend actually only on A' ' 
and not on A) : 

f»(»o»»-|^-tr(^ + ^jg)) = -*»8,(^) + 
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where fx L0 = ^ - iA<® and fi R0 = £t(H^ + *A( ))E. 

ill it! 

The introduction of the above shorthand notations /i^o and /i^o is motivated by the fact that they can be identified 
as chemical potentials but for the left and right chiral densities (they are thermodynamic conjugate variables of pl 
and pft). Nevertheless, due to the fact that we are considering the SU(2) case and no more £7(1), new features appear, 
that is we have now some couplings with the derivatives of these chemical potentials /irq and /iLo- 

On the right-hand side of the equation, the entropy production is, as usual, a bilinear form between "thermodynamic 
forces" and dissipative fluxes. Again, we can make linear combinations with those thermodynamic forces and introduce 
transport coefficients. Our choice for the hydrodynamic velocity imposes also that u^j^ = and u^j R = and we 
have to make the entropy production positive. These prescriptions allow us to eliminate some coefficients and to 
obtain some constraints on the remaining others. We get : 

K = (a,id^ + [CskA (f^£t)^ Ti (21) 

= or - u^n (n dv + « £ ,,.a, (^) + KRA . dv (^)) (22) 

f L = - (K Lji d v (g) + [K LL ] id d v + O^j ) n (23) 

f R = Or - ( ( ^ + [KLilli A + [«H*kA f ^ ))n , 24 ! 



( 3 ^ - uV) ({C2 - ^v)dxu x + Ctjd x (^i(Zd x tf)X\ + 77 ((g» x - u»u x )d x u» + {g Xv - u x u v )d x u») (25) 



where [Q] means that Q is a 3x3 matrix. Due to the Onsager reciprocity principle, all the matrices except [klr] 
are symmetric, 4£i.j = £4^ and there are actually 39 independant coefficients. If we represent the quadratic form 
of the entropy production by a 12 x 12 matrix, we can easily show that all coefficients appearing in the diagonal 
should be positive and that there exists inequalities between the 39 coefficients : as for the superfluid, all the principal 
minors have to be positive. We see that dissipative equations for SU(2)l x SU(2)r symmetry implies new transport 
coefficients and couplings between baryonic, left-handed and right-handed currents. We also see that spontaneous 
breaking of the chiral symmetry implies the existence of matricial transport coefficient for, by instance, the bulk 

viscosities £1 and C3 : the thermodynamic force d\ ^^i(S9 A E^) : ,^ which is the equivalent of d\{V 2 d x 4>) has now 

three components. These components are nevertheless all multiplied by the same factor 4f- which is related to the 
amplitude of the order parameter of the chiral symmetry. Of course, near the transition phase this amplitude is not 
frozen as in our case and can fluctuate. The dynamics of these fluctuations is such that the amplitude of the order 
parameter has to be incorporated explicitely in the approach since then it becomes an hydrodynamic variable \l2\ . 

Finally let us remark that, since it is known that dissipation can affect the observables (see |l3j for the modification 
of the temperature profile used to describe heavy ion collisions), it will be of fundamental importance to determine 
quantitatively the influence of these new couplings. 



CONCLUSIONS AND OUTLOOK 



We have treated in this paper the dissipation in the hydrodynamic regime for relativistic systems with spontaneously 
broken symmetry U(l) and SU(2). For the superfluid, we recovered the non-relativistic limit of 0. For the SU{2) 
case, we introduced new transport coefficients associated to the right and left charges and to the Goldstone modes. 
Since we know that transport coefficients can affect quantitatively observables, it would be interesting to express them 
with Kubo-type relations and then compute them explicitely from this microscopic approach. This work is under 
study |3| ■ An other extension of this paper is to determine the new relaxation times with the effect of the symmetry 
breaking in order to know if this implies some noticeable modifications to the conclusions drawn in 

3] 

concerning the 
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typical equilibration time of heavy ions collisions. It is well known that the way we introduced the dissipative effects 
in that paper leads formally to instabilities and can not allow to reach the underlying physics of the relaxation times. 
In order to solve that problem one has to introduce them explicitely. One convenient way to do this is to use the 
14-moments of Grad as done in yl. 
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